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SUMMARY 

The deposition rate of boric oxide from a stream containing submicron 
particles was found experimentally to fit a molecular diffusion equation 
as follows: 

Particle diffusion rate (particles/(sq ft) (sec)) 

= 18.08 (velocity of particles treated as ideal gas molecules (ft/sec)) 
times (particle volume concentration in free stream (particles/cu ft)) 

The particle-velocity expression, taken from the kinetic theory of 
gases, assumes that the particle acts like a molecule of an ideal gas of 
molecular weight equal to the product of the particle mass times Avogadro's 
number. The particle volume concentration is calculated from knowledge 
of the fuel, fuel-air ratio, temperature, and pressure. 

This equation enables the calculation of boric-oxide deposition rates 
to within ±15 percent on surfaces where conventional flow patterns exist. 
The assumptions and conditions limiting the use of the equation are dis- 
cussed. The range of the experimental program covered typical turbojet 
primary-combustor operating conditions. 

A discussion of the effect of boric oxide deposits on jet engine 
performance is included. 


INTRODUCTION 

The deposition of boric oxide B 2 O 3 on surfaces in the turbine region 

of jet engines operating on boron- containing fuels results in engine- 
performance deterioration (refs. 1 and 2). A better understanding of the 
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mechanism of "boric oxide deposition in this region could result in reduc- 
tion of these losses, and thus would hasten the time when "boron hydrides 
would "be used operationally as fuels for turbojet engines. For this 
reason, an investigation of the deposition characteristics of "boric oxide 
was performed at the NACA Lewis laboratory. 

Reference 3 suggests that the first step in solving this problem is 
the determination of growth rates of liquid-boric oxide particles. Data 
on these growth rates are presented in reference 4. The particle diameter 
for residence times and conditions existing in typical turbojet engine 
primary combustors was found to be between 1.0X10"^ and 1.5X10 ^ centi- 
meters. It was also suggested in reference 4 that boric oxide deposition 
is a diffusional process because of the extremely small particle sizes. 

By using this theory as a starting point, the original analysis was modi- 
fied; and a theoretical equation, expressing deposition rate as a function 
of exhaust -gas stream parameters, was derived. 

The results of both the theoretical and experimental work are re- 
ported herein, along with a discussion of the effect of deposition on 
engine performance. 

This report is divided into three parts. The first (THEORY) gives 
in detail two possible mechanisms for the deposition of near molecular- 
size (1000- to 2000-angstrom units) particles on a surface and develops 
theoretical equations for each. The second part (EXPERIMENTAL APPARATUS 
AND PROCEDURE, RESULTS AND DATA CORRELATION, AND CALCULATION METHOD) 
presents the results, shows how the necessary empirical constants were 
obtained to compare properly the results with theory, and applies the 
results to a method of calculating boric oxide deposition rates on engine 
surfaces. The third part (DISCUSSION AND APPLICATION OF RESULTS, AND 
CONCLUDING REMARKS) discusses the applicability of the calculation method 
and the effect of deposits on engine performance. 

An example of the use of the calculation method appears in appendix 

B. 


THEORY 

General Discussion of Problem 

The theoretical analysis is based on the results of the particle- 
size determinations of reference 4; that is, the particles are small 
enough to assume that their motion is essentially molecular.. From this 
assumption, a net particle flux from regions of relatively high concentra- 
tion to regions of low concentration might occur by two processes: mo- 

lecular diffusion, turbulent diffusion, or a combination of both. (Dif- 
fusion caused by Brownian motion also might be considered, but it is 
negligible for the conditions of interest.) 
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In general, when a stream is turbulent, diffusion processes occur 
within the boundaries of the stream at a much faster rate than that pre- 
dicted by molecular diffusion considerations. However, the present dis- 
cussion involves particles which are small enough that their motion is 
essentially molecular, but which are quite large when compared with typi- 
cal gas molecules such as oxygen or nitrogen. If the scale of turbulence 
in the stream approached the diameter of the particles, an expression 
describing the deposition of these particles on surfaces might be expected 
to contain terms corresponding to both molecular and turbulent diffusion. 
Therefore, the two diffusion mechanisms will be treated separately, and 
the experimental data will be used to determine which one, or combination, 
actually is controlling the deposition process. 


Turbulent Diffusion 

In order to obtain a turbulent diffusion rate, the well-known 
Reynold’s analogy for the relation between the transfer of momentum and 
the transfer of mass in a boundary layer on a surface will be used. This 
relation can be stated as follows: 

Rate at which mass is transferred to surface 
Rate at which mass flows in free stream 

. Rate of momentum loss in boundary layer due to friction 
Total momentum of free stream flowing over surface 

which can be written 


n t = 

hqUo P q Uq 


( 1 ) 


(A list of the symbols used in this report is given in appendix A.) 
In equation (l), the only factors that are not immediately obtained from 
a knowledge of free-stream conditions are the turbulent diffusion rate 
Nrp and the shear stress at the surface t w . 


Reference 5 gives an expression for x w , provided that the shape of 
the velocity profile in the boundary layer is assumed. A good approxi- 
mation in most applications is to assume that 


U = U 0 (y/8) 1//7 = UqII 1 / 7 


( 2 ) 


where 


n = (y/&) 


( 3 ) 
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By using equation (2) and the definition of t w in reference 5, it can 
be shown that 


T 

w 
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g etc 



(l - q 1 / 7 )!) 1 / 7 
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(4) 


Normally, the Boundary -layer thickness can he expressed as follows: 
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(5) 


Then 
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= k(l-a)(^ x-a = 


k(l-a) 

Re a 


Substitution of equation (6) into equation (4) gives 


U W 


7 P0 U 0 k(l-a) 
= 72 S R e a 


( 6 ) 


(7) 


Substitution of equation (7) into equation (l) gives an expression for 
the turbulent diffusion rate N T in terms of known stream parameters: 


N t = 0.09722 n 0 U Q k hl|) 

Re 


( 8 ) 


Molecular Diffusion 

Loeb (ref. 6) shows that the flux of molecules N passing through a 
unit area in unit time can be expressed by 


B - i 5 n 


(9) 


Now, if one mean free path away l (in the y direction) the concen- 
tration is n-z(^), the net flux in the y direction is 


N 


net 


- s - i [n 


7(dn/dy)] 


( 10 ) 
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which reduces to 


N 


net 



an 

dy 



or, hy defining the diffusion coefficient D as 



(ID 


( 12 ) 


then 


N m ~ D dj 


(13) 


which expresses the molecular diffusion rate in terms of a known diffusion 
coefficient and a concentration gradient. 


Discussion of Theoretical Equations 

Generally, the diffusion of a molecular species to or from a turbu- 
lent stream is governed by the turbulent diffusion rate. However, in the 
application discussed in this report, the liquid particles, although 
assumed to be governed by the laws of molecular motion, are rather large 
as compared with normal molecular sizes. The scale of turbulence decreases 
rapidly near a surface; therefore, the scale of turbulence might approach 
the particle size as the surface is approached. This could mean that, 
while turbulent diffusion might aid in bringing the particle near the sur- 
face, molecular diffusion actually would be the driving force across a 
very thin film on the surface, which might be called a laminar sublayer 
with respect to the particle (see ref. 7). Because of the physical prop- 
erties of the particles, it is assumed that every particle colliding with 
the surface remains at the surface. Therefore the particle concentration 
approaches zero at the wall. Based on this discussion, the particle con- 
centration profile might be illustrated by the following sketch: 

Velocity profile Particle concentration profile 
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EXPERIMENTAL APPARATUS AND PROCEDURE 
General Apparatus 

The experimental program consisted of two phases: the first, to 

determine the "boundary -layer characteristics of the test surface; and the 
second, to determine the deposit characteristics of the test surface. 

The apparatus used in each phase was the same except for the 
instrument at ion . 

A photograph and a schematic diagram of the test apparatus are shown 
in figure 1. Combustion air from the laboratory central air-supply system 
entered the apparatus, was measured with an ASME standard orifice, and 
was throttled through a remotely operated throttle valve. The air then 
entered the inlet plenum and went through a series of wire and cheese 
cloth screens that removed dirt and insured a uniform velocity profile. 

The combustor (just downstream of the inlet plenum) is a standard jet- 
engine tubular combustor. The fuel injector is air atomizing and is the 
same as reported in reference 8. From the combustor, the exhaust products 
entered a calming section where the velocity was about 60 feet per second 
or less. In this section, any large oxide particles sheared from the 
combustor walls settled out; thus the diameter of all the particles enter- 
ing the test section would be between 1.0*10~5 and 2.0*10“5 centimeters. 

The test section (fig. l) consisted of a duct of constant width (4 in.) 
and varying height (3.70 in. at the inlet and 18.20 in. at the exit). A 
wedge, 4 inches wide and 12 inches long with upper and lower surfaces 
parallel to the top and bottom of the test section, was used as the test 
surface. The wedge was positioned in the foreward part of the test sec- 
tion, and a rectangular spacer filled up the remaining length. The upper 
surface of the wedge was the test surface. Two ducts carried the gas 
stream from the test section into the exit plenum, where the gases were 
quenched with water sprays and were discharged into the building exhaust 
system. 

The fuel systems for pentaborane and JP-4 fuel were the same as 
reported in reference 9. The fuel system used in the trimethyl borate - 
methyl alcohol azeotrope (TMlB) tests consisted of a 2000-gallon storage 
tank, a pump, remotely operated throttle and shutoff valves, and provisions 
for purging the fuel lines. A diagram of the TMB fuel system is shown in 
figure 2. The fuel flow in each system was measured with a rotating-vane- 
type flowmeter. 


Instrumentat ion 

The following table summarizes the instrumentation that was used in 
both the boundary -layer and deposition tests: 
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Station 

Station description 

Number of tem- 
perature points 

Number of 
pressure points 

1 

Combustor inlet (air) 

3 

1 (static) 


Combustor inlet (fuel) 

1 

0 

2 

3 

Combustor exit 

Upper , lower (or both) 

7 (av.) 

1 (static) 


test-section exit ducts 

5 

1 (static), av. 
of 5 (total) 


During the "boundary -layer tests, the test surface was instrumented 
with surface thermocouples and surface static -pres sure taps (12 of each) 
placed 1 inch apart along the test surface. 


Boundary -Layer Test Procedure 

The Boundary -layer thickness and velocity profile were measured with 
a total-pressure prohe. The details of the prohe are shown in figures 3 
and 4. The pressure impulse was converted to a voltage by a pressure 
transducer, and the output was recorded on two instruments: (l) the Y- 

axis of an X-Y self-balancing recording potentiometer, which made visual 
monitoring of the boundary -layer profile possible while a run was in prog- 
ress, and (2) a high-speed, self-balancing recording potentiometer that 
was used because of its greater accuracy. The probe was moved by a re- 
motely operated actuator equipped with a variable precision resistor that 
indicated the probe position. The output from the resistor was recorded 
on the X-axis of the X-Y recorder. 

When the probe, electrically insulated from the rest of the apparatus, 
touched the test surface, it closed a circuit that actuated the pen on the 
X-Y recorder, which in turn indicated the position of the test surface in 
relation to the actuator travel. The total actuator travel was l/2 inch, 
and the accuracy of the position indicator was h^D. 003 inch. The pressure 
transducer and the position indicator were calibrated before and after 
each test . 


Deposition Test Procedure 

In making the deposition tests, inlet airflow and temperature were 
established, and the combustor was ignited on JP-4 fuel. The combustor 
outlet temperature was adjusted to the test conditions. When temperature 
equilibrium on the test surface was reached, either the TMB or pentaborane 
fuel flow was initiated, and the test fuel flow was adjusted. When the 
test fuel had entered the combustor, the JP-4 fuel was turned off, and 
the conditions were held constant for some predetermined time at which 
the run was terminated. 
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The upper wedge surface (test surface) was washed to remove the 
boric oxide deposit, and the resultant boric acid solution was chemically 
analyzed for the total quantity of boric oxide. A mean value of the depo- 
sition was taken because, after a very short running time, some flow of 
the oxide occurred along the surface and made accurate determinations of 
local deposition rates impossible. In tests where oxide was blown back 
into the separation zone (discussed in a later section) , the wedge and 
the separated surfaces were analyzed separately. 

The residence time was calculated by using flow-continuity relation- 
ships and the dimensions of sections between the combustor and the test 
section. From the residence time, a particle size in the test section 
was calculated by using the method of reference 4. 


Several preliminary tests were conducted in order to determine the 
test-section velocity and temperature profiles, the temperature drop be- 
tween the test section and the exit ducts, and the time required for the 
test section to come to temperature equilibrium. The results of these 
tests showed that about 8 minutes of warm-up time was required. The local 
velocity across the test section varied less than 3 percent from the 
average for the range investigated. Typical test-section temperature 
profiles are shown in figure 5. 


Fifty-nine boundary -layer thickness measurements were made over the 
following range of conditions: 


The coordinate system used in the boundary -layer studies is illustrated 
in figure 6. 

The results are presented in table I and figure 7. A line was fitted 
to the data, which resulted in the following equation: 


Preliminary Tests 


RESULTS AND DATA CORRELATION 


Boundary -Layer Studies 


Free-stream velocity, ft/sec 

Temperature, °R 

Reynolds number (xUq/Vq) . . 


1.50x105 to 4.71X10 6 


44 to 1983 
498 to 1618 



( 14 ) 
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From the boundary- layer data, the probable error of any data point refer- 
enced to equation (14) was calculated to be ±5.58 percent. As seen in 
figure 7, the data fall somewhat below the classical flat-plate line. 

This result was not unexpected, since reference 7 shows this trend for 
surfaces with a slightly positive pressure gradient (as on a wedge) . 


Deposition Studies 

Eighteen deposition tests were made, thirteen with TMB and five with 
pentaborane. The data are summarized in table II. Photographs of the 
wedge and test section after typical trimethyl borate (run 9) and penta- 
borane (run 15) tests are shown in figure 8. As seen in figure 8(a), the 
deposit on the wedge surface was quite uniform up to the point where the 
wedge met the spacer block. Flow separation occurred at the point where 
the wedge and spacer block met and resulted in two effects: (l) in the 
region of flow separation, the deposition rate was greatly reduced; and 
(2) boric oxide, pushed from the wedge into the separated region, remained 
there and continued to build up (fig. 8b). 


The first step in the data correlation was the introduction of the 
constants a and k from the experimental boundary -layer data into equa- 
tion (8) (turbulent diffusion equation), which gives: 


Nm 

— = 0.0105 
n 0 


U, 


0 


0.22 

Re 


(15) 


If turbulent diffusion is the controlling mechanism, a plot of (N t t /n n ) 

0 2 2 

against (U 0 /Re’ ) would be expected to result in a straight-line log-log 

plot, although the slope might be other than unity. This plot was made, 
and no line could be drawn through the data. A second plot of N tot 
0 22 

against (U 0 n 0 (Re‘ ) was made, and the following equation was taken from 

data (fig. 9) : 


"tot - 5 - 51xl ° 


6/ u Q n Q 


,0.514 


vRe 


0.22 


(16) 


This equation indicates that turbulent diffusion is not the control- 
ling mechanism for, if it were, N-tot would be directly proportional to 
n 0 . The terms involved in a molecular diffusion process were also examined 
by determining an effective diffusion coefficient Dg-pp and a distance 

which this diffusion coefficient was controlling. This distance will be 
referred to as the diffusion boundary -layer thickness and will be denoted 
by 6 d . In determining D 0ff and 5 D , the following assumption was made: 
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The concentration of particles at y = 5^ is equal to the free -stream 
concentration uq. This assumption, which is valid only in cases where 
uniform flow (no separation) is present, is justified "because the flow of 
particles through the diffusion boundary layer, rather than the turbulent 
diffusion through the turbulent boundary layer, will control the total 
flux to the wall. 

If an expression similar to equation (ll) is used to describe the 
diffusion process, an explicit solution for both D e ff and cannot 

be obtained by using the experimental data available. Therefore an as- 
sumption regarding one of the two unknowns must be made, and it was assumed 
that the diffusion boundary layer is identical in thickness to the mean free 
path in the binary mixture (ref. 6, p. 99). In other words. 


By using this assumption and the experimental data, the effective diffu- 
sion coefficient was calculated from the following equation: 

D = M ~tot Z 0 

eff 

The effective diffusion coefficient was compared with the theoretical 
molecular diffusion coefficient (using the mean free path of the binary 
mixture and assuming the particle velocity to be equal to the theoretical 
molecular velocity) and was found to vary linearly with it, as indicated 
by the following expression: 




D eff = 54 * 23 D m ( 19 ) 

which is illustrated in figure 10. The absolute arithmetic average devia- 
tion of all the experimental data from equation (19) was ±12.5 percent. 

The implication is that molecular diffusion across a thin layer in the 
boric oxide - air mixture controls the boric oxide deposition rate. 

The combination and rearrangement of equations (12), (18), and (19) 
give the following experimentally determined expression for the diffusion 
of boric oxide particles to a surface: 


H tot ' 18 - 08 Vo 


( 20 ) 
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CALCULATION METHOD 


Derivation 


The following discussion illustrates how equation (20) can be used 
to calculate boric oxide deposition rates. The assumptions used in the 
derivation of the equations necessary for the solution of the deposition 
rate will be consistent with the assumptions used in reference 4, which 
are: 

(1) The liquid particle may be treated as a large ideal gas molecule 

(2) The particles are spherical. 

(3) The particle concentration is constant across the main flow 
stream. 

(4) The free-stream conditions are such that essentially all of the 
boric oxide is in the liquid phase. 

In order to calculate deposition rates on surfaces , the parameters 
n Q and £p must be known. The following discussion shows how the par am 
eters are found and how they can be used to calculate deposition rates. 

Loeb (ref. 6) shows that the average velocity of an ideal gas mole- 
cule may be expressed by 



(21) 


For the case of particles acting as ideal gas molecules, 


m = W p N Av 


( 22 ) 


where 



(23) 


Therefore 



( 24 ) 


12 


MCA RM E57F18 


Before solving equation (24), dp must "be calculated. Reference 4 shows 
that 

dp = (o.2056 liia *VST W ‘ (25) 

In this equation, only v 1;L /v (the ratio of liquid volume to gas 
volume in the system, which is a function of fuel type, cp, P, and T) is 
unknown . 

The number of pound-moles of gas phase per pound of boric oxide X 
formed during combustion (assuming no dissociation or BgOj vapor), along 
with the operating temperature, static pressure, and boric oxide liquid 
density (p 1;L = f(T) only), can be combined to yield v liq / v g in the 

following manner: X = f(fuel type, equivalence ratio) only. Let 

Y = P/p llq RT (26) 

Then 

liii = 1 (27) 

v g -X 

which equates ( v iiq/ v g) to f ue l type, cp, T, and P. If v liq/ v g •'- s 
known, the particle diameter may be found, and then 



which gives all the information necessary for a deposition-rate 
calculation. 


General Method 

(Each equation is dimensionally correct and will yield a correct, 
result as long as consistent units are used.) The steps for calculating 
the deposition rates are as follows: 

(l) Start with the fuel type and fuel-air ratio to calculate the 
equivalence ratio cp ; 


Actual fuel-air ratio 
Stoichiometric fuel-air ratio 


<P 


(29) 
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(2) By knowing the fuel type and <P, get X from figure 11. 

(3) By knowing the temperature and static pressure, get Y from 
figure 12. 

(4) Calculate ( v iiq/ v g ) from 

(v.^ /v ) = Y/X (dimensionless) (27) 

(5) Find residence time of particles from formation to point where 
deposition rate is desired using flow continuity relationships and areas 
of the system. 

(6) Calculate particle diameter from 

s = (°- 2056 ^ t V ^)°' 4x10 ' 4 (25) 

(r = d yq e . - £l — T = °K; P-, . = gm/cm 3 ; t = sec; and d = cm] 

\ (mole) (°K) xiq y ) 


(7) Calculate n^ from 


n<- 


v g 


= v -m ( -4 


Ttdr 


(28) 


(8) By knowing T and d p , calculate ? p from 


£ p = 2.84X1CT 12 


RT 


f p, . d 
M liq p 


(24) 


(It is more convenient to use c.g.s. units in equation (25) and to 
convert to feet per second by dividing by 30.48 cm/ ft) 

(9) Calculate N tot from 


»tot = 18 ‘ 08 Vo 

(10) Calculate the deposition rate r from 

r = N tot (g - rtd p) 


( 20 ) 


(30) 
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where r is given in feet per second if N^^ is stated as particles 
per square foot second and dp is given in feet. 

An example of the calculation of boric oxide deposition on a typical 
turbine stator blade is shown in appendix B. 


DISCUSSION AND APPLICATION OF RESULTS 
Deposition Tests 

The results of the deposition tests indicate that the preceding cal- 
culation method may be used to calculate boric oxide deposition rates on 
engine surfaces where conventional turbulent flow patterns exist (turbine 
stator blades and annular walls , tailpipes, exhaust nozzles, and ramjet 
engine combustors) . The analysis is not valid in some regions that exist 
in jet engines: for example, surfaces where flow separation occurs, and 

surfaces on which a film of clean air is introduced (through pores or 
louvers) . On surfaces where flow separation occurs, the deposition rate 
has been observed to be considerably less than on surfaces where uniform 
flow is present (see fig. 8 and discussion in previous section) . The data 
of reference 8 indicate that deposits are greatly reduced on surfaces 
where air filming is employed (in ref. 8, a turbojet combustor), making 
calculation of deposition rates on these surfaces unnecessary. 

The deposition rate alone, however, does not indicate completely the 
problem caused by the presence of liquid combustion products. If deposi- 
tion occurred (with no run-off) at the rate of 4*10“4 inches per second 
(a typical figure for a pent aborane -fueled turbojet engine), a thickness 
of 1.44 inches would be present on the engine surfaces at the end of 1 
hour of operation. In actual practice, deposits do not become this thick. 
The liquid oxide flows along surfaces because of the shear stress exerted 
on it by the flowing gas stream. At points where abrupt changes in surface 
profile occur, the flowing liquid film is sheared from the surface into 
the gas stream. After a period of time, the deposition rate equals the 
run-off rate; an equilibrium condition is then reached, and the film thick- 
ness will be a constant (as long as operating conditions remain constant) . 
Although no data are presently available on the equilibrium film thick- 
ness, it is logical to assume that the thickness would be a function of 
either the shear stress acting on the liquid film or the kinetic force 
acting on waves formed in the film, similar to one of the following models: 
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At the equilibrium condition (where the run-off rate is equal to the 
deposition rate), the deposition rate is important because it can be used 
in determining the amount of oxide running off turbine stators or being 
thrown from the tip of turbine rotors by centrifugal force. 


Turbojet Engine Performance 

Boron hydride fuels apparently show the greatest advantage in the 
primary combustor of the turbojet engine. Unfortunately, this application 
comprises the major problems arising from boric oxide deposition, and 
these difficulties occur mainly with the turbine (refs. 1 and 2). Shifts 
in turbine performance (which in turn can cause a fall-off in over-all 
engine performance) are shown in at least four ways : 

(1) Formation of deposits on turbine stator blades reduces the flow 
area and thus requires an increased turbine- inlet pressure to maintain 
constant engine mass flow. The increase in turbine-inlet pressure (or 
compressor -discharge pressure) necessitates an increase in turbine work 
at constant mass flow. 

(2) Any increase in turbine pressure ratio (increased turbine-inlet 
pressure or decreased turbine-discharge pressure) at constant engine mass 
flow results in a closer approach to the condition of limiting loading of 
the turbine rotor blades. As this condition is approached, the turbine- 
discharge Mach number increases, which results in high tailpipe total- 
pressure losses (accentuated by boric oxide deposits on tailpipe surfaces) 
and shock losses downstream of the throat section of the rotor blade. 

(3) The presence of boric oxide deposits on turbine surfaces reduces 
the aerodynamic efficiency of the turbine because of increased friction 
losses . 

(4) Changes in turbine operating conditions (such as increased pres- 
sure ratio) can cause losses in other engine components. For example, a 
higher turbine-discharge Mach number M^ can increase the tailpipe 
total -pressure loss due to friction (this friction loss is proportional 

to M?,). The effect of these losses on the performance of an actual en- 
td 

gine is illustrated in the following discussion. 

In the full-scale engine test of reference 1, a J65 engine (2-stage 
turbine) was run on pentaborane for 11 minutes. The engine was warmed up 
on 62-octane gasoline and it showed a specific fuel consumption (sfc, 
pounds of fuel per pound of net thrust per hour) of 1.365. Immediately 
after introduction of the pentaborane, the specific fuel consumption 
dropped to 0.942 (an improvement of 31 percent) but increased with run 
time until, at the end of the test (11 min) , a value of 1.240 was reached. 


16 


NACA RM E57F18 


However , since engine conditions did not remain constant during the run, 
the specific fuel consumption after 11 minutes of operation cannot he 
compared directly with that of gasoline at the beginning of the test. 

This decay in engine performance was due mainly to an increase in 
tailpipe total-pressure loss which was caused by an increase in the 
turbine -discharge Mach number. This increase in turbine -discharge Mach 
number could have been caused by one or both of the following mechanisms : 

(1) During the test, the turbine -stator area was calculated to have 
decreased about 5 percent (from 0.757 to 0.720 sq ft) because of boric 
oxide deposits. Consequently, to maintain constant engine mass flow, the 
turbine-inlet pressure had to be increased; this required an increase in 
compressor work and a corresponding increase in turbine work (which in- 
creases the turbine pressure ratio) . The net effect was an increase in 
turbine -discharge Mach number M^a, which increased the friction pressure 
drop in the upstream portion of the tailpipe (this loss is proportional 
t° M^ d ). 

(2) After the test was completed, a buildup of boric oxide was ob- 
served near the trailing edge of the suction side of the second- stage 
stator blades. This buildup is attributed to flow separation that proba- 
bly occurs near the trailing edge of these blades . The force of the gas 
stream acting on the oxide film on the forward part of the blade pushes 
oxide into the separation zone, where there is very little force present 
to remove it. This buildup of oxide pushes the point of flow separation 
upstream along the blade. The effect is a decrease in the tangential 
velocity vector at the stator exit, which requires a greater pressure 
ratio for the same turbine work output. To keep turbine work constant, 
the absolute velocity vector at the stator exit must be increased by in- 
creasing the turbine pressure ratio; this increases the turbine-exit Mach 
number (fig. 13) . 

This oxide buildup in zones of flow separation also occurred in the 
test section used during the tests reported herein. Figure 8(b) shows 
that a large accumulation of oxide occurred on the spacer-block surface 
just downstream of the wedge. Flow separation occurred at the intersec- 
tion of the wedge and spacer block, and the oxide deposited on the wedge 
was pushed into the separated region, where the force on the oxide to 
produce further flow was very slight (fig. 14) . 


CONCLUDING REMARKS 

Experimental data indicating that molecular diffusion across a thin 
layer is the controlling factor in the deposition of boric oxide are 
presented. A diffusion equation containing an experimentally determined 
constant is given, which enables the calculation of boric oxide deposit 
rates to within ±15 percent on surfaces where conventional flow patterns 
exist . 
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Knowledge of the behavior of liquid exhaust products is not complete. 
Additional work on run-off rates is required for the calculation of equi- 
librium thickness. Combined data on deposition rates and equilibrium 
thickness would greatly aid in the understanding of deposit patterns and 
subsequent engine performance losses. Such information also would be of 
value in alleviating the deposit problem. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 5, 1957 
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APPENDIX A 
SYMBOLS 

constant defined by eq. ( 5 ) 

diffusion coefficient, sq ft/sec 

particle diameter, cm or ft 

acceleration due to gravity, 32.17 ft/sec^ 

constant defined by eq. (5) 

mean free path, ft 

Mach number 

molecular weight 

particle flux due to diffusion, particles/(sq ft) (sec) 
Avogadro T s number, 6.023x10^ molecules/mole 
particle volume concentration, particles/cu ft 
static pressure, atm 

universal gas constant, units consistent with rest of equation 

Reynolds number , (xU 0 / v ) 

deposition rate, ft/sec 

absolute temperature, °R or °K 

time, sec 

velocity, ft/sec 

local velocity, ft/sec 

volume, cu ft 

weight, gm 

number of lb -moles of gaseous combustion products/lb of B2O3 
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x,y,z distance along coordinate axis, ft 

Y liquid volume concentration parameter, Y = p/p RT 

liq 

8 full boundary-layer thickness, ft 

bp effective diffusion boundary- layer thickness, ft 

T) defined by eq. (3) 

v kinematic viscosity, sq ft/sec 

£ molecular velocity, ft/sec or cm/sec 

p density, Ib/cu ft 

t w wall shear stress, Ib/sq ft 

cp equivalence ratio 

Subscripts : 

eff effective 

g gaseous phase 

liq liquid phase 

m molecular 

net net 

p particle 

T turbulent 

td turbine discharge 

tot total 

0 free-stream conditions 
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APPENDIX B 

CALCULATION METHOD USED TO DETERMINE DEPOSITION RATE ON 
TYPICAL TURBINE STATOR BLADE 

The stator -"blade shape chosen for this calculation was taken from 
reference 10. The Mach-number distribution for this particular blade 
shape is shown in figure 15. 

The calculation was made for the following conditions: 


Fuel 

Equivalence ratio 

Turbine-inlet static temperature, °R 

Turbine-inlet static pressure, atm 

Combustor-exit Mach number (turbine inlet) . . 

Ratio of specific heats 

Stator -outlet mean Mach number • 

Stator Mach number distribution (from ref. 10) 
Particle diameter, cm 


Pentaborane 

.• . . . 0.2 

. . . 2700 

. ... 1.5 
. . . 0.325 
. . . 1.30 

. . . 0.915 
. . Fig. 15 
. 1.2xl0' 4 5 6 


The results of the calculation are summarized in table III and on 
figure 16. 
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TABLE I. - BOUNDARY- LAYER DATA 


C\D 

CO 


Run 

Temper- 

ature, 

T, 

°R 

Dis- 

tance 

along 

axis, 

x, 

ft 

Reynolds 

number. 

Re 

Bound ary- 
layer 
thick- 
ness, 5, 
ft 

6/x 

(a) 

Run 

Temper- 

ature, 

T, 

°R 

Dis- 

tance 

along 

axis, 

x, 

ft 

Reynolds 

number. 

Re 

Boundary- 
layer 
thick- 
ness, 6, 
ft 

6/x 

(a) 

1 

498 

0.833 

4 . 707X10® 

0.00451 

0.00542 

31 

523 

0.583 

3.45X10 6 

0.00328 

0.00563 

2 





4.423 

.00392 

.00470 

32 

520 



3.288 

.00333 

.00571 

3 





3.405 

.00372 

.00447 

33 

515 



3.018 

.00305 

.00523 

4 





2.381 

.00458 

.00550 

34 

505 



3.726 

.00290 

.00497 

5 





5.272 

.00400 

.00480 

35 

505 



2.55 

.00300 

.00514 

6 

701 



1.493 

.00516 

.00620 

36 

504 



1.72 

.00543 

.00931 

7 

745 



2.076 

.00512 

.00615 

37 

531 

.333 

.150 

.00327 

.00982 

8 

750 



2.521 

.00425 

.00510 

38 





.300 

.00292 

.00878 

9 

770 



2.834 

.00416 

.00500 

39 





.532 

.00223 

.00699 

10 

773 



3.232 

.00425 

.00510 

40 





.574 

.00262 

.00787 

11 

765 



3.428 

.00401 

.00481 

41 

1591 

.500 

.8143 

.00340 

.00679 

12 

1249 



1.441 

.00541 

.00650 

42 

1587 



.8157 

.00352 

.00703 

13 

1548 



1.186 

.00487 

.00585 

43 

1546 



.9953 

.00350 

.00701 

14 

1367 



1.505 

.00650 

.00780 

44 

1607 



1.0992 

.00377 

.00754 

15 

1560 


r 

1.371 

.00517 

.00621 

45 

1618 



.9511 

.00327 

.00654 

16 

508 

0.417 

2.072X10 6 

0.00215 

.00515 

46 

791 

0.500 

0.09822X10 6 

0.00530 

0.01060 

17 

505 



1.112 

.00301 

.00723 

47 

763 



.1417 

.00483 

.00966 

18 

501 



1.696 

.00251 

.00603 

48 

742 



.1851 

.00486 

.00973 

19 

501 



2.040 

.00230 

.00552 

49 

517 

.833 

.4709 

.00526 

.00632 

20 

500 



2.959 

.00218 

.00524 

50 





.2679 

.00724 

.00869 

21 

498 



2.882 

.00205 

.00493 

51 





1.165 

.00501 

.00602 

22 

706 



1.631 

.00245 

.00588 

52 





.7204 

.00584 

.00701 

23 

700 



1.498 

.00251 

.00601 

53 





.5399 

.00800 

.00961 

24 

717 



1.225 

.00282 

.00676 

54 

642 



1.911 

.00419 

.00503 

25 

720 



.859 

.00296 

.00711 

55 





1.595 

.00499 

.00599 

26 

1210 



.218 

.00425 

.01020 

56 





1.321 

.00518 

.00622 

27 

1416 



.298 

.00341 

.00818 

57 





.7871 

.00583 

.00700 

28 

1349 



.281 

.00355 

.00852 

58 





.5446 

.00636 

.00763 

29 

1346 



.251 

.00375 

.00900 

59 





.4519 

.00659 

.00791 

30 

556 

\ 

f 

1.75 

.00382 

.00917 










a 6/x, ratio of boundary-layer thickness to distance from leading edge. 
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TABLE II. - SUMMARY OF DEPOSITION TESTS 


Run 

Fuel 

Temper- 

ature, 

OR 

Pres- 

sure, 

atm 

Test 

section 

veloc- 

ity, 

ft/sec 

Particle 

concentra- 

tion, 

particles 

Measured 

deposition 

rate, 

particles 

Mean free 
path of 
particles, 
ft 

Average 
velocity 
of par- 
ticles, 
ft/sec 

Theoretical 
molecular- 
diffusion 
coefficient, 
sq ft/sec 

Effective 
diffusion 
coefficient, 
sq ft/sec 

cu ft 

sq ft-sec 

1 

a TMB 

1912 

1.41 

901 

4 . 22 X 10 13 

2 . 54 X 10 14 

1 . 87 X 10- 15 

3 . 58 X 10-* 1 

2 . 23 X 10- 15 

1 . 13 X 10" 14 

2 


1878 

1.11 

758 

3.62 

2.48 

2.56 

3.68 

3.14 

1.75 

3 


1740 

1.38 

813 

4.05 

1.98 

1.84 

3.48 

2.13 

0.899 

4 


1716 

1.07 

658 

3.34 

2.04 

2.52 

3.57 

3 .00 

1.54 

5 


1631 

1.36 

855 

4.58 

2.50 

2.23 

3.73 

2.78 

1.22 

6 


1655 

1.13 

926 

5.17 

3.76 

3.78 

4.33 

5.45 

2.75 

7 


1560 

1.37 

746 

4.11 

2.49 

2.04 

3.58 

2.43 

1.24 

8 


1572 

1.07 

588 

3.20 

1.88 

2.60 

3.58 

3.10 

1.53 

9 


1410 

1.46 

685 

4.07 

2.48 

1.90 

3.52 

2.23 

1.16 

10 


1490 

1.13 

510 

2.83 

1.67 

2.08 

3.30 

2.28 

1.23 

11 


1360 

1.17 

461 

2.74 

1.61 

1.95 

3.21 

2.09 

1.13 

12 


1287 

1.02 

552 

3.70 

3.00 

3.70 

3.96 

4.88 

3.00 

13 


1320 

1.02 

658 

4.43 

3.75 

4.54 

4.34 

6.56 

3.84 

14 

b B 5 H g 

1716 

1.00 

690 

2 . 99 xl 0 13 

1 . 46 xl 0 14 

0 . 882 xl 0" 15 

2 . 21 xl 0 _1 

0 . 649 xl 0’ 15 

0 . 432 xlO -14 

15 


1884 

1.04 

1123 

4.99 

2.66 

1.64 

2.96 

1.61 

.870 

16 


1530 

1.01 

798 

4.17 

3.30 

1.79 

2.96 

1.76 

1.41 

17 


1749 

1.01 

1004 

4.85 

2.74 

2.05 

3.19 

2.18 

1.16 

18 


1841 

1.02 

984 

4.33 

1.94 

1.45 

2.78 

1.34 

0.652 


a TMB, Trimethyl borate - methyl alcohol azeotrope. 
b B5Hg , Pentaborane . 
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TABLE III. - CALCULATION OF TURBINE STATIC BLADE DEPOSITION RATES 


CO 


Sta- 

tion 

Distance 

from 

leading 

edge, 

ft 

Pound- 
moles of 
gaseous 
combustion 
products/ 
lb of 

Tem- 

pera- 

ture, 

T, 

OR 

Static 

pres- 

sure, 

p 

. 9 

atm 

Liquid 

volume 

concen- 

tration 

parameter, 

v 

Ratio of 
liquid to 
gaseous 
phase 
volume, 

V liq/ V S 

Particle 

concen- 

tration, 

n Q , 

particles 

Par- 

ticle 

veloc- 

ity, 

ft/sec 

Total par- 
ticle flux 
due to 
diffusion, 

N tot > 
particles 

Deposition 
rate , 
ft/sec 





1 

cu ft 





B2O3, x 





(sq ft/sec) 







Suctioi 

n surface 






0.01667 

.03334 

.05000 

.06668 

.08333 

.1000 

.1167 

.1333 

.1500 


0.01667 

.03334 

.05000 

.06668 

.08333 

.1000 

.1167 

.1333 

.1500 


0.82 


2615 

1.31 

0 . 84 x 10 

2557 

1.20 

.71 

2412 

.908 

.47 

2332 

.784 

.439 

2321 

.767 

.424 

2352 

.812 

.453 

2376 

.851 

.466 

2412 

.908 

.47 

2431 

.939 

.519 


.-5 


1 . 02 X 10 
.866 
.573 
.535 
.517 

.552 
.568 
.573 
.633 


-5 


3 . 18 X 10 14 

2.54 

1.78 

1.67 

1.61 

1.72 

1.77 

1.78 
1.97 


0.771 
.760 
.737 
. 723 
.721 

.729 
.731 
.737 
. 740 


4 . 43 X 10 15 

3.49 

2.37 

2.19 

2.10 

2.26 

2.34 

2.37 

2.64 


1 . 42 X 10 
1.12 
.761 
. 704 
.673 

.726 

.750 

.761 

.846 


-4 


Pressure surface 


0.82 


2711 

1.51 

0 . 75 x 10 ' 

2716 

1.52 

.75 

2711 

1.51 

.75 

2695 

1.48 

.72 

2670 

1.43 

.71 

2637 

1.36 

.67 

2596 

1.26 

.63 

2524 

1.13 

.61 

2459 

.990 

.45 


^5 


0 . 915 X 10" 5 

.915 

.915 

.878 

.866 

.817 

.768 

.744 

.549 


2 . 85 X 10 14 

2.85 

2.85 

2.74 

2.70 

2.54 

2.39 

2.32 

1.71 


0.783 

.783 

.783 

.781 

.779 

.772 

.766 

.753 

.743 


4 . 04 X 10 15 

4.04 

4.04 

3.87 

3.80 

3.55 

3.31 

3.17 

2.30 


1 . 29 X 10- 4 

1.29 

1.29 

1.24 

1.22 

1.14 

1.09 

1.02 

.739 


Leading edge 


0 

0 

0.82 

2700 

1.50 

0 . 72 X 10" 5 

0 . 878 X 10 -5 

2 . 74 X 10 14 

0.781 

3 . 87 X 10 15 

1 . 24 X 10 -4 
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To exhaust 



Test surface 
Test section 


Combustor 


Calming section 


Throttle valve 


Orifice 


Combustion air 


. Inlet plenum Fuel injector 
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A-A 


Exit plenum 


Water spray 
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Figure 1. - Test apparatus. 
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Figure 3. - Detail of total-pressure tube used in boundary- 
layer studies. (All dimensions are in inches.) 
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Figure 4. - Wedge and test section showing "boundary- layer 
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Figure 5. - Typical test-section temperature 
profiles . 
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Figure 6. - Coordinate system used in boundary- layer and deposition tests 
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Ratio of boundary-layer thickness to distance 
from leading edge, 6/x 
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(a) A trimethyl borate fuel. 
Figure 8. - Oxide deposits on wedge. 
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Figure 8. - Concluded. Oxide deposits on wedge 
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Figure 9. - Preliminary correlation of deposition data. 
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Effective diffusion coefficient, D e ff, sq ft/sec 
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Theoretical molecular diffusion coefficient, D m , sq ft/sec 


Figure 10. - Comparison of experimentally determined diffusion coefficient with theory. 
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(a) Trimethyl "borate - methyl alcohol azeotrope fuel. 

Figure 11. - Pound-moles of gas phase per pound of "boric oxide 
in combustion products plotted against equivalence ratio. 
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Pound-moles of gas phase per pound of boric oxide in combustion 

products, X 
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Equivalence ratio, <P 


(b) Pentaborane fuel. 


(c) Diborane fuel. 


Figure 11. - Continued. Pound-moles of gas phase per pound of boric oxide in combustion products plotted 
against equivalence ratio. 




NACA EM E57F18 


6 


CD 


& 

0 
•H 
•P 
03 

3 

1 

o 

o 


'S 


o 


u 

o 

P 

ft 

o x 

Tj «S 

a w 
3 p 
o o 
ft d 

?H O 

0 u 

ft ft 
<u 

03 

3 

J3 

ft 

CQ 

a) 

bO 

ft 

0 

03 

<D 

1 

i 

Ti 

d 

3 

o 

ft 


• 1 


.2 


(d) Ethyl decaborane fuel. 
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(e) Propyl pentaborane fuel. 




Figure 11. - Concluded. Pound-moles of gas phase per pound of boric oxide in combustion products 
plotted against equivalence ratio. 
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Figure 12. - Boric oxide liquid-volume concentration parameter as function of temperature and pressure. 
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0-0^ Change in direction angle of 
absolute velocity vector due 
to deposit 

U Absolute velocity vector 
without deposit 
Absolute velocity vector 
with deposit 

U Axial velocity vector 
without deposit 

U z £ Axial velocity vector 
9 with deposit 

U^-U Change in absolute velocity 
vector due to deposit, nec- 
essary to hold Ug constant 

Ug Tangential velocity vector 

(constant for constant turbine 
work) 



Figure 13. - Effect of flow separation on boric oxide buildup on trailing edge of turbine stator 
blades, with resulting effect on absolute velocity vector. 
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Figure 14. - Boric oxide buildup in region of flow separation on test surface. 
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Free -stream Mach number adjacent to 
boundary layer 
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Figure 15. - Stator Mach number distribution. Exit Mach number, 
0.915. 
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